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ABSTRACT 


This  paper  analyzes  the  regularity  of  plaine  eleisticity  problems  with 
piecewise  analytic  data.  The  results  are  given  in  the  fr2unes  of  the  countably 
weighted  Sobolev  spaces  introduced  by  BabuSka  emd  Guo  [SIAM  J.  Math.  Anal., 

Vol  19  (1988),  pp.  172-203,  Vol.  20  (1989),  pp.  763-781].  These  results  are 
of  major  significance  for  the  practical  design  and  theoretical  analysis  of  the 
p  auid  h-p  version  of  the  finite  element  solutions  of  the  elasticity 
problems. 

Key  Words.  Elasticity  problem  with  piecewise  euialytic  data,  countably  normed 
space,  weighted  Sobolev  space,  singularity. 


AMS(MOS)  Subject  Classifications.  35^5,  35D10,  35G15,  35J05,  73C02 


§1.  INTRODUCTION 


The  reguleu'ity  of  the  boundary  value  problems  for  elliptic  equations  was 
addressed  in  meLny  papers.  The  regularity  of  the  solution,  when  the  domain  is 
not  smooth,  was  addressed  in  [18,19]  zuid  subsequently  in  many  papers  (see  e.g. 
[13,20,21,22,26]).  The  auialysis  was  made  for  the  scalar  equation  (e.g. 

Laplace  equation)  and  for  general  elliptic  systems,  e.g.  theory  of  elasticity. 
Although  the  results  are  essentially  similar,  the  proofs  are  in  various 
details  different.  The  regularity  of  the  solution  is  typically  described  in 
terms  of  Sobolev  spaces  and  various  decompositions  of  the  solution  in  the 
singular  and  regular  parts.  These  results  are  important  not  only  for 
theoretical  reasons  but  also  for  the  design  and  euialysis  of  the  numerical 
treatment  of  these  problems. 

Finite  element  method  is  today  the  mostly  used  method  for  solving  the 
elliptic  equations.  There  au'e  three  versions  of  the  finite  element  method, 
the  h,  p  and  h-p  finite  element  method  (see  e.g.  [7]).  The  following 
questions  arise  in  the  connection  with  the  FEM: 

a)  What  is  a  typical  class  of  problems  in  applications,  e.g.  structural 
mechanics, 

b)  What  reguleu'ity  results  are  the  most  useful  for  the  design  and 
ainalysis  of  the  finite  element  method,  especially  its  h-p  version. 

The  answer  to  the  first  question  is  that  in  the  structural  mechsinics  we 
deal  with  the  problems  with  piecewise  analytic  data  (boundzu'y  of  the  domain, 
coefficients  of  the  equation  and  boundary  conditions,  etc.).  Description  by 
the  standard  Sobolev  of  finite  order  neglects  features  important  for  the 
numerical  solution. 

The  answer  to  the  second  question  is  that  we  need  such  a  description 
which  allows  to  construct  the  numerical  method  which  is  maximally  effective. 
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In  the  case  of  the  p-version  of  the  finite  element  the  user  provides  the 
mesh,  emd  the  degree  is  adaptively  selected  (see  e.g. ,  Applied  Structure 
(Rasna  Corp.  CA,  USA)  or  STRIPE  (Aeronautical  Resezu'ch  Inst,  of  Sweden)).  The 
proper  regularity  description  leads  to  the  rules  which  the  user  has  to  follow 
(see  e.g.  [27]).  In  the  case  of  the  h-p  version  the  proper  form  of  the 
regularity  results  allow  to  construct  the  meshes  and  degrees  such  that  the 
convergence  of  the  method  is  exponential,  see  e.g.  [4,5,6,14]. 

We  have  found  that  the  proper  form  of  the  regularity  results  which  serves 
best  the  goals  of  the  numerical  auialysis  is  the  description  of  the  regularity 
in  the  frame  of  countably  normed  (weighted)  Sobolev  spaces.  Although  the 
spaces  zu'e  close  to  the  Gevrey  spaces  (see  [9,11])  they  have  some  special 
different  features,  and  the  form  of  the  statements  au'e  directly  related  to  the 
effective  use  in  the  FEM.  We  have  introduced  these  spaces  in  [2,3,4]  and 
shown  their  slgnificzint  applications  for  the  finite  element  analysis.  We  have 
proven  that  the  h-p  version  with  an  exponential  rate  of  convergence  csui  be 
constructed. 

In  [2,3,4]  we  addressed  the  problem  of  the  scalar  equation.  Here  we 
address  the  elasticity  problem.  Although  the  bzisic  results  and  the  analysis 
methodology  are  similar  here  and  for  the  scalar  problem,  there  au'e  various 
essential  details  which  need  significant  modification  of  the  auialysis. 

Because  of  the  major  importance  of  the  elaisticity  in  application,  the  detailed 
theoretical  results  presented  here  au'e  essential.  They  do  not  follow  directly 
or  easily  from  the  available  results  in  the  literature. 

In  Section  2  we  introduce  the  weighted  Sobolev  spaces  and  countably 
normed  spaice.  The  vau'iational  solution  of  the  elaisticity  problem  with  data 
given  in  weighted  spaces  is  addressed  in  Section  3.  Section  4  gives  a 
complete  auialysis  of  the  elasticity  problem  in  an  infinite  angulau'  domain. 
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In  Section  5  we  prove  the  regularity  theorems  for  the  elasticity  problem  in 
polygonal  domains  in  terms  of  the  countably  weighted  Sobolev  spaces. 

The  methodology  and  techniques  used  for  the  elasticity  problem  can  be 
used  for  general  elliptic  systems  of  equations. 


3 


§2.  PRELIMINARIES 


The  notations  and  definition  of  spaces  in  [14]  are  the  ones  most  often 

used  in  this  paper.  Let  Q  denote  a  polygon  with  vertices  and  open 

edges  connecting  A^  and  1  <  i  <  =  A^^)  shown  in  Fig.  2.  1 

Let  2)  auid  M  be  subsets  of  M  =  {1,2, . . .  ,M} ,  7)  =M\M,  ^ 

ie2) 


U  r  r  =  an  =  r°  U 

ie// 

Let  3  =  (3j . 3^) 

be 

M-tuple  of 

real  numbers, 

0  <  3^  <  1,  1  S  i  <  M,  and 

where 

k  is  an 

integer  8uid 

r^  is  the  distance  between 

X  and  vertex  A^. 

Fig.  2. 1  Polygonal  Domain 

By  H*'(n)  we  denote  the  usual  Sobolev  spaces  and  by  the 

P 

weighted  Sobolev  spaces  for  integer  k,  t,  k  >  ^  >  0  with  the  norm 

k 


llu||\  -  =  llull^.  ,  + 

iC’^n)  r'^(n) 

p 


lal^e 


(0) 


(if  1^0,  the  term 


a  »  (aj.a2).  a  0 


2  (X 

llulP-  -  drops  out)  where  Du 

r"^(n) 

integer,  i  =  1,2,  We  shall  write 


a'«'u 

“l  “2 
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For  convenience,  we  use  the  spaces  ^(Q)  instead  of  H^’^(Q)  if  Q 

P  P 

is  finite  or  infinite  sector,  with  the  norm 


(Q)  r‘^(Q)  L'^CQ) 


a'“'u 


(if  i  =  0,  the  term  ||u||  .  drops  out)  where  2)  u  =  u  =  - 

a,  a_  a,  a 


2  -  1=0  2 
r  0  5r  30 


(r,0)  is  the  polar  coordinate,  a  is  the  same  ais  above. 

The  spaces  W^(Q)  introduced  by  Kondrat6v  (see  [16,  17,  18]),  are  used 
P 

in  this  paper  as  well  in  polar  coordinates,  and 


,  ,,2  V*  r  2(j3-k+ai)  ,^a  ,2  , 

lull  .  =  >  r  ^  12)  ul  rdrd0. 

w:(Q) 


Let  D  =  {(t,0)|,  -<»<t<«,  0<8<ta>,  suid  for  integer  k  i  0  auid 


real  h  >  0  we  define 


h{J(D)  =  (u|  Y  f  e^^'*'|D“u|^dTd0  =  ||u||\  <  col- 

We  shall  write  H^(D)  =  L^(D). 

I 

The  countably  normed  spaces  B_(n),  t  =  0,1,2,  are  defined  by 

P 

8^(0)  =  (ulu  €  H^’^w),  ll♦^^,^_^D“u||^  <  Cd^^’^k-i)!, 

for  |a|  =  k  =  ^, <+l, . . . ,  with  C,d  i  1  independent  of  k> , 

I 

and  ®^(S)  is  defined  in  jxjlar  coordinates  on  finite  sector  S  =  ((r,0)| 
0<r<5,  O<0<w}  by 

S^(S)  =  (ulu  €  i^'hs).  ||r^~^*“»2)“u||^^  S  Cd*''^k-«)!. 


for  |a|  =  k  =  with  C,d  2:  1,  independent  of  k} . 

We  shall  write  Bf({2, C,d)  or  sf(S,C,d)  if  we  emphasize  the  constants  C 
P  p 

amd  d. 
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and 

12)^u|^  =  Y  r‘^“2  |2)“u|2 

lai^k 

The  theorem  on  the  equivalence  of  weighted  Sobolev  spaces  defined  in 
Cartesi2ui  coordinates  and  polar  coordinates  will  be  used  later.  We  quote  the 
following  theorem  from  (2). 


Theorem  2. 1  (cf.  Theorem  2.1  of  (2]).  Let  S  be  a  finite  sector  =  {(r,0)| 
0<r<5,  O<0<u>,  *13  =  0  <  3  <  1,  then  for  0  <  ^  <  2  integer, 

u€H^’^(S).  (resp.  B^(S))  if  and  only  if  u  e  '^(S)  (resp.  sf(S)). 

P  P  P  P 

We  denote  the  vector  and  vector  space  by  bold  face.  For  example,  u  = 
(Uj.Ug)^.  H*'(n)  =  H^(n)  X  B*(n)  =  B^(n)  x  B^(n),  etc.  The  theorem 

quoted  above  holds  for  vector  spaces  ^(S)  amd  (resp.  B^Cn) 

P  P  p 

auid  sfcn))  as  well. 

P 
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§3.  VARIATIONAL  SOLUTION  OF  ELASTICITY  PROBLEM  WITH 
DATA  GIVEN  IN  WEIGHTED  SPACES 


(3. la) 


Consider  the  linear  plain  strain  elasticity  problem  on  polygon  Q 

■  -mAu,  -  dl  •  Ig)  .  f, 

)  -  a  3  a  in  £1 

-MU,  -  d-g  *  ig]  - 


(3. lb) 


1  0  ,u0, 

“'rO  *  '  “  ° 

T  =  <r-nl  ,  =  g'  •  g'i 


where  u  =  (Uj,U2)  is  displacement,  and  <r 
tensor  given  by 


L®’21.0'22j 


is  the  stress 


(3.2) 


=  2n 

dui  ^ 

3U2- 

11 

dXi 

'^Xi 

3X2- 

=  24 

3U2' 

22 

3X2 

(^Xi 

3X2* 

r3ui  ^  ^ 

2-] 

®^12 

^  ®‘21 

L  "  ^ 

'^X2  3X] 

iJ  • 

A  and  n  are  Laun^  coefficients,  n  is  unit  outward  normal  to  dtJ,  the  force 
f  =  (fj,f2)^.  the  boundary  condition  =  (gj.g^)^,  =  (G^.G^)^,  etc. 

Using  the  differential  operator  L  and  the  boundary  operator  B  we 
write  (3.1)  as 


(3.  r ) 


Lu  =  f 

Bulj,  *  Ig°,g^] 


in  a 


Theorem  3. 1.  Let  f  €  L„(n)  =  H°’°(n),  g^  e  ^^^"^(F^) ,  ^  =  0.  1,  and 

P  P  p 

ir^l  *  0,  then  problem  (3.1)  has  a  unique  solution  u  €  H^(n)  (in  the  weak 
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sense ) ,  and 


(3.3) 


Proof.  First  we  assume  that  =  0.  The  bllineau'  form  on  Hq(Q)  x  HQ(n) 
,3.4,  1^)  ,  .  (1^  .  g)  [1^  *  1^] 


+  M 


to  >  £^l]dxdy 

[axt  axzj  [axi 


where  H^(iJ)  =  (ulu  €  H^(Q),  ul^o  =  0^- 

Owing  to  Lemma  2.9  auid  2.11  of  [2]  we  have 

I  f  g^.vdsl  <  C  llg^l  llvtl 

Jfi  „y2,l/2(j.lj  ^1 

P 

From  Lemma  2. 10  of  (2]  we  have  f  €  (H^(CI))',  auid 


(0) 


,  V  V  €  H^(n) 


i|  f*vdxi  s  c'  iifiij^  iivii^^j  ,  V  V  €  r(n), 


'Q 


“0^”^  H"(n) 


Hence 


F(v)  =  f  f»vdx  +  f  g^* 

Jn  •'r‘ 


V  ds 


is  a  llnesu'  and  continuous  functional  on  H^(n),  and 

1,. 


IIFII 


i  C 


(H"(n))' 


[liriL 


.]• 


The  variational  problem  is  to  seek  u  e  Hq((2)  such  that 
B(u,  v)  =  F(v)  for  auiy  v  e  H^Cn). 

There  is  a  well-known  Korn  inequality,  i.e. ,  for  u  €  11^(0).  IF^I  *  0, 

">'"'^,0,  ‘  ^  ^  '“"h> 


(n)  ’ 


is 
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where  >  0  (i  =  1,2)  are  indepjendent  of  u  (see  [17]).  Applying  the 
Laoc-Milgram  theorem,  we  have  existence  and  uniqueness  of  the  weaJc  solution  of 

(3.1),  and  (3.3)  holds.  If  *  0,  let  w  =  u  -  G*^,  then  w  e  HQ(n) 
satisfies 

B(w,  v)  =  J  f‘vdx  +  J  ^  ^g^  -  l^j'vds  -  B(G^,v),  V  v  e  ^^(£1) 

=  F(v). 

Obviously  F  e  (H^(n))'.  Applying  the  result  above  we  obtain  the  existence 
and  uniqueness  of  the  solution  of  (3.1)  in  general  (in  the  wesiJc  sense). 


§4.  AUXILIARY  PROBLEMS  ON  INFINITE  SECTOR 


In  order  to  analyze  the  behavior  of  the  solution  at  the  corners  of  Q,  we 
consider  the  elasticity  equation  (3.1)  on  am  infinite  sector  Q  = 

{(r,0)|O  <  r  <  «,  0  <  9  <  (j> .  We  prefer  to  write  the  equation  in  polar 

f  = 


coordinates  (see  e.g.  (251).  Let  u  =  (u  ,u„)  and  <r 

r  0 


j^rr '  *^r0'j 


(i  etc.,  and  let  A  -  ("fl 

r  0  l^-sin0,cos0j 

T  - 

=  AcrA  .  u  satisfies 


-  -  -I  I  - 

Then  u  =  Au,  f  =  Af,  g  =  Ag  ,  <r 


(4.  1) 


-d[4u,  -  ^  -  (A  .  M)  i  »  jX  «  1  , 


in  Q, 


amd  <r  is  given  by 


rr 


2“  3F^  •  *  r  *  F  ae^] 
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(4.2) 


c.  3u-  u  >  fSu  u  ,  3u.-, 

-ee  >  *  f"]  *  •  h-/] 


d\i  au^ 


<r  =  <r 
re  0r 


1  _£  +  _0  -  _0l 

r  30  3r  r  J 


One  of  three  kinds  of  boundary  conditions  may  be  imposed  on  3Q,  namely 


(4.3a)  u  1  =  =  G^l ,  ^ 

0=0 ,  u)  *  0=0 ,  u 


(Dirichlet  condition) 


=  (<r^Q.<rgQ)'^  =  =  G^lg^o.o)  (Neumann  condition) 


“  <0=0  =  ^  =  ^'0=0  ’  ^®'r0’‘^00^^>0=u,  =  «i  =  ‘^^'9=0,  condition) 


—t  —I  —t  —t  -I 

Hereafter,  g  denotes  the  different  vectors  g-  =  G  (r,0)  and  g  =  G  (r.u) 

in  the  trace  sense  at  0  =  0,«.  Let  L(2))  and  B(2))  be  corresponding 

matrix-differential  operator  and  boundaury  operator.  We  rewrite  (4.1)  and 

(4.3)  as 


(4.4) 


(L(2)).  B(2))]G  =  (f,  i°,  i^]. 


By  the  change  of  vairiable  x  =  £n  p  .  we  convert  the  equation  (4.1)  and 
boundary  condition  (4.3)  into  the  problem  L(u)  =  f  on  a  strip  domain  D  = 
{(t,0)  I  -<»<T<«,  O<0<  u},  l.e.  , 


(4.5) 


2  2 
^  u  ^  >  3  u 

-(2n+A)  ~  M  ^  +  (p  +  A) 


3^u  3u 

3t^  +  (3p  +  A)  ^  =  f^ 


in  D 


2~  ~  2~ 

3  u  3u  3  u 

aia?  tI  -  "  rr 

OX 


b\. 


-  (2n  +  A) 


-  +  liu  =  f 

0  ^0  0 


with  one  of  the  following  boundairy  conditions: 
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(4.6a)  u  I  -  =  =  G*^Iq  n  •  (Dirichlet  condition) 

W— U  f  U)  o— U » (*) 

(4.6b)  T  =  (®’.,a. ®’qq)^  =  n  •  (NeumaJin  condition) 

To  oo  o— U,{i> 


(4.6c) 


~  I  _  ~0  =0, 

“  ^0=0  ~  *0  "  ®  1 


0=0,  ^  ~  ^'^T0 


=  g^.  = 


90 


0— C()  U)  0*w 

(Mixed  condition) 


where 


(4.7) 


u  =  =  (u^(e~’*^,0), 

.0))’'',  £  =  0,1,  f  =  (?.^.?q)^ 


u  (e''^.0))'^,  =  (G^,G^)^  =  e"^^(G^e"'^.0). 

o  To  P 

=  e'^'^(f^(e"'^.0),  fg(e“^,0))’^’ 

/^u  >  ^  3u 

'ee  *  [ae-  *■  "tJ  *  ^  *  “t  * 


Further,  by  Fourier  transformation  we  obtain  a  system  of  ordinau-y 
equations 


(4.8) 


2“^  A 

-n  — 5^  +  (2h+X)(1+tj^)u  +  (i7}(M+A)  +  (3fi+X))  33-  =  f 

d0'^  T  d0  T 


du 

-((3ji+A)  -  iT)(/i+A))  -  (2m+A) 


2  - 

—  .  M(U,  )u. 


with  one  of  the  following  boundary  conditions: 


(4.ga) 

u 

*0=0, w 

■0  AOi 
-  *  =  G  1, 

(4.gb) 

A 

T| 

A  A 

=  (<r  <r 

'0=0, w 

T0'  00 

(4. gc) 

ul 

0®O 

"0  -0, 

*0  ■  ®  'e- 

0=0, w  ’ 

,T. 


(Dirichlet  condition) 


''1  ''1 

=  g  =  G  I-  _  ,  (Neumann  condition) 

o— U » (j) 


‘1  “1, 


-  |0=W 

(Mixed  condition) 


where  u  »  (u^,u  )^  »  (?(G),  ^(u-))  ?(u),  etc., 

T  o  V  T  o 
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?(u)  =  — ^  r  e  ^’*^u(T,0)dT,  T)  =  ^  +  ih,  ^  h  >  0, 

'^-00 

and 

A  A 

“  /<iu  ^  ^  ^  du 

(4.10)  -  (Itj+DUqJ  .  o-gg  =  (2n+X)  ^  +  (2M+A-iAT,)u^. 

By  L(D  , 7))  amd  B(D  , tj)  we  denote  the  matrix-differential  opjerator  and 

tj  & 

boundary  operator.  Furthermore,  we  write  the  system  and  the  boundary 
condition  as 

(4.11)  U(v)u  =  (L(Dg.T}).  B(Dg.T,)]u  =  [f.G°,G^]  (or  [f .  ] .  ) 

The  operator  IKt})  =  [L(Dg,T)),  B(Dg,T))]  depends  polynomially  on  the 
complex  parsuneter  tj.  By  the  argument  used  in  [14]  for  all  tj,  with  the 
exception  of  certain  isolated  points,  H(tj)  realizes  an  isomorphism: 

H^(I)  »  L^(I)  X  H^(I)  X  H^(I)  (or  L^(I)  x  x  C^).  Consequently,  the 

inverse  operator  X(v)  ®  IKij)  ^  is  an  operator- valued  meromorphic  function  of 

TJ  with  poles  of  finite  multiplicity.  These  poles  are  the  eigenvalues  of 

IKtj)  (see  [10,16,18,19,20,21,22]).  For  each  pole  tj  of  3?(tj),  the 

homogeneous  problem  of  (4.11)  has  at  leaist  one  non-trivial  solution 

corresponding  (eigenvector  function)  in  1^(1).  The  transcendental  equations 

which  the  eigenvalues  satisfy  have  been  derived  in  several  different  ways  zind 

can  be  seen  in  the  literatures  of  continuum  mechanics  euid  mathematics  (see 

e.g.  [10,  24,  26,  28]).  The  typical  approach  is  to  consider  a  bihairmonic 

equation  in  Q  insteaul  of  the  elasticity  equations.  Since  we  adopt  the 
-  T 

displawrement  u  =  polau*  coordinates  in  (4.8)  and  (4.9)  the 

A  A 

coefficients  of  the  operator  L(Dg,Tj)  and  B(Dg,Tj)  80*0  constaints.  Hence  we 
are  able  to  derive  the  transcendental  equations  directly  from  the  homogeneous 
equation  (4.8)  and  the  boundairy  condition  (4.9)  on  displeicement  and  traction. 
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Since  these  equations  are  not  of  the  main  interest  in  this  paper  we  will 
present  them  in  a  lemma  below.  For  completeness  we  include  the  proof  in  the 
Appendix. 


Lemma  4.1.  Let  tj  =  iz  be  an  eigenvalue  of  U(i}) ,  and  u  be  the  Poisson 
ratio.  Then 

(i)  for  the  Dirichlet  problem  (4.8)  and  (4.9a).  z  satisfies 

2 

(4.  12)  sin^zu  =  ^3^4^! 

(ii)  for  the  Neumann  problem  (4.8)  and  (4.9b),  z  satisfies 

(4. 13)  sin^zu  =  z^sin^u  , 


specially  ti  =  0  is  an  eigenvalue  with  multiplicity  of  2,  the  corresponding 

—  T  —  T 

eigenvectors  are  =  (cos  0.  -sin  0)  and  =  (sin  0,  cos  0)  , 


(iii)  for  the  mixed  problem  (4.8)  and  (4.9c),  z  satisfies 
(4. 14) 


,  2  4(l-u)^  z^  ,  2 

sin  zu  «  -7=— 3 - =—7—  sin  u 


3-4u  3-4u 

and  7}  =  ±i  (i.e.,  z  =  ±1)  are  the  eigenvalues  with  multiplicity  of  1  if 
(1  +  — )  cos  2u  +  1  =0. 


From  the  equation  (4.11)  ~  (4.13)  it  is  easily  seen  that  zeroes  of  these 
equations  are  symmetric  with  respect  to  the  origin  euid  the  real  aixis  in 
complex  pl2Lne.  Hence  the  eigenvalues  of  TKt))  eu'e  located  in  the  complex 
plane  symmetrically  with  respect  to  the  origin  and  the  imaginary  axis.  By  T 

V 

we  denote  the  eigenvalues,  and  let  tc^  be  a  positive  number  such  that 

(4.15)  K.  =  min  Imi)  -  min  |  ImT)|  . 

T)€T^  TJ€T^ 

V  V 

Im  T)>0  Im  71*0 

Next  we  prove  the  Agranovich  &  Vishik  conditions  I  and  II  which  are 
substantial  to  the  key  inequality  (4.17)  (see  [1]).  These  two  conditions  were 
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used  implicitly  or  explicitly  in  many  papers,  e.g. ,  [261  for  elasticity 
problems  and  in  [21,22]  for  general  elliptic  systems. 

Let  Du  =  i  ^  auid  A^CD, t))  be  the  principal  part  of  the  operator 

A  A 

L(Dg,T)).  We  write  Aq(D,7})  in  matrix  form 


Aq(D,  tj) 


Ti^(2|i+X)  + 
(ji+A)T}D 


(fi+A)T}D 

M7»^  +  (2fi+A)D^ 


Lemma  4.2.  (Condition  I)  For  ^  (real),  t)  e  Z .  =  {tjI  larg  tjI  <  p,  or 

Pi  1 

larg  T)  -  111  <  with  any  e  (0,ir/2)  amd  ItjI  +  |^I  *  0,  det  (Aq(^,t})) 

*  0.  Furthermore,  the  equation  det  (Aq(<,t)))  =0  in  ^  has  equal  numbers 

of  roots  in  upper  zuid  lower  half-planes  for  t)  6  Z,  and  t)  *  0. 

9x 


Proof.  It  is  easy  to  see  that 

det  (Aq(^,7)))  =  p(2m+A)(7>^+C^)^  *  0 

for  <  €  R^,  TJ  €  Z^^  with  any  P^  €  (0,jr/2),  and  1^1  +  ItjI  *  0.  Also  it  is 

A 

seen  that  ^  =  ±  irj  are  the  roots  of  the  equation  det  (Aq(^,t)))  =0  in  C 
(complex).  Hence  the  equation  has  2  roots  in  upper  and  lower  plane,  respec¬ 
tively  if  0  *  TJ  e  Z;.  .  o 

A  ^ 

Let  Bg(D,Tj)  be  the  principal  part  of  the  boundau'y  operator  B(Dg,Tj) 
defined  by  (4.10),  and  0=0  (resp.  «),  I  =  (0,o8)  (resp.  I_  =  (-<»,  u)). 

U  9  0 

Then  we  have  the  following  lemma. 


Lemma  4. 3.  (Condition  II)  For  any  p  e  (0,ii/2),  if  tj  *  0  2uid  tj  e  Z, 

X  Pi 

=  (tjI  larg  11  -  Til  <  pj^  or  larg  tjI  <  p^},  the  equation  on  the  half  line 


(4.  16) 


Ag(D,,)W  -0.  8  < 
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has  a  unique  stable  solution  W  such  that  |W|  — >0  as  0— (resp. 
0  — »  -») . 


Proof.  We  will  prove  the  lenuna  for  0^  =  0  and  =  (O.oo).  The  proof  for 

u  ©0 

0-  =  w  is  similar  with  what  follows.  For  the  homogeneous  equation  (4.16) 


the  solution  W  must  have  the  form  e 


with  b  satisfying  the 


equation 


det  (AQdb.Tj))  =  M(2p+X)  (b^-t)^)^  =  0. 

Hence  b  =  ±tj  are  the  roots  with  multiplicity  of  2.  For  t)  g  S  and  rj  * 

<Pi 

0,  Re  b  =  Re  T)  *  0.  Let  a  =  -sgn  (Ret)),  cr  = 


fw  1  fi 

u  _  ^  ai?0  ^  c 

We  .1 


0  +  fxjor/ri 


is  the  stable  solution  if  Cj^  and  can  be  uniquely  determined  by  any 

*  2 

given  boundary  condition  h  €  C  . 

A  A 

For  the  Dirichlet  condition  W|_  _  =  h,  c,  and  c.,  satisfy 

1 0— a  1  2 


1  a<r/'i)l  fc 


=  h  . 


Then  obviously  c^  and  can  be  uniquely  solved  from  the  equation  above. 

For  the  Neumann  condition  we  have  by  (4.10) 


Bq(D,tj)W  = 


M  ^  -  iMTjWg 


t-iAi,W^  +  (2p+X) 


Then  the  Neumann  boundary  condition  B(D, i))W  =  h  leads  to 


'2iin  M(l+<r)  fcj 

i2M7|  i(2M+A(l-<r))  °2 


=  h  . 


IS 


Since 


det 


■  ’2^7) 

pd+o-) 

■ 

=  i2p7)  (X+p)  (l-<r) 

=  -i4p^7)  *  0 

.Ii2p7) 

i(2p+A(l-<r)) 

• 

j  euid 

Cg  can  be  uniquely  determined  by 

h.  Therefore  (4.16) 

hEis  unique  stable  solution  for  tj  €  Z .  and  t)  *  0  smd 

9\ 


-boa 


— >  as  0  — >1 


with  Cq  >  0  and  =  -|Re7)|  . 

After  verifying  the  conditions  I  and  II  we  have  the  following  theorem. 


Theorem  4.1.  Suppose  there  is  no  pole  of  on  the  line  Imi)  =  b,  then 

A  A  A  ^  A  4 

the  solution  u  *  31(7)) [f,G  ,G  ]  of  the  problem  (4.8)  and  (4.9)  satisfies 

(4.17)  IlGlI^  +  |7)|^||u||^  i  C  fllfll^  +  y  IIG^ll^ 

iTd)  L'^d)  ^  L'^d)  ir’^d) 

t=o,  1 


Proof.  Due  to  Lemmais  4.1  and  4.2  the  conditions  I  and  II  eu'e  satisfied  on  an 
angle  =  {i\:  larg  t)|  <  or  |arg  t)  -  nl  <  with  €  (0,7r/2).  By 

Theorem  6.1  of  [1],  (4.17)  holds  with  C  independent  of  t)  and  u  if  t)  e  Zj 
and  |7)|  >  7)q,  where  7)^  is  some  positive  real  constant.  The  line  Im  t)  =  b 
is  contained  in  Z^  except  a  finite  segment  for  which  IRe  7)|  <  |b|  cosec 
Hence  (4.17)  holds  for  v  on  the  line  Im  7)  =  b  with  |7)|  >  7)^.  For  those 
7)  on  the  line  with  |7)|  S  7)^,  31(7))  is  analytic.  Hence  for  f  €  L^(I),  g^  e 


c^. 


the  solution  u  exists  in 


H^d), 


suid 


Hull 


+  |7)|^|G||%  s  c  (1+I7)i*)(llfll%  +  y  diil^+lg' 

ITd)  L'^d)  ^  L^d)  ^  " 


')) 


^=0,  1 
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ir(i) 


I 

e=o,  1 


"t  2 

IT  ‘(I) 


^  2-  1^1  ('«o'  ^’«u' 

1=0.  1 

where  aaid  are  some  constants  independent  of  u  and  t).  Thus  we 

have  proved  (4.17)  for  t)  on  the  whole  line  Im  t)  =  b  on  which  K(t))  has  no 
pole.  □ 


Lemma  4. 4.  Let  7  be  a  strip  =  {Tjl-h  <  Im  <  h>,  0  <  h  <  k^,  then  3?(7))  has 
no  pole  in  J  for  the  Dirichlet  euid  mixed  problems,  and  the  origin  is  the 
only  pole  of  JKtj)  in  7  for  the  Neumann  problem. 


Proof.  Due  to  the  definition  of  and  Theorem  6.1  of  (ll  may  have 

poles  on  a  finite  segment  of  the  real  axis.  We  shall  show  that  there  is  no 
pole  of  JKtj)  on  the  real  axis  for  the  Dirichlet  and  mixed  problems  amd  the 
origin  is  the  only  pole  of  JiCij)  on  the  real  aixis  for  the  Neumauin  problem. 

(i)  Suppose  that  v  ~  ^  (real)  is  the  pole  of  3?(7j)  for  the  Dirichlet 
problem,  ^  satisfies 

.  .  sin  w  - 

‘  3^  «  ■ 

But  Lemma  A.  1  shows  that  for  C  =  0  u  =  kir,  there  is  no  non-trivial 

solution  of  (4.8)  and  (4.9a).  Hence  we  may  eissume  that  u  *  kn.  Let  f(^)  = 
shu)^  -  |sin  C4)l^,  then  f' (?)  »  uch?u  -  |sin  a>|  S  w  -  Isin  ul  >  0.  Hence  f(?) 
>  0  for  0  <  ?  <  OB.  Simllau'ly,  we  can  show  that  f'(?)  <  0  for  -oo  <  ^  <  0. 
Therefore  |sh?(i>|  >  Isinwl  |?|  >  |sin  ul  |?|  for  real  ?  *  0  and 

V  6  (0,1/2),  which  implies  that  the  zero  is  the  only  real  number  which 
satisfies  (4.12).  Since  zero  is  not  an  eigenvalue,  31(ti)  has  no  pole  on  the 
real  line  for  the  problem  with  the  Dirichlet  boundary  condition. 
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The  proof  for  the  mixed  boundary  condition  is  similar. 

(ii)  Suppose  tj  =  C  (real)  is  a  pole  of  for  the  Neumann  problem.  Then 

^  satisfies 

sh^o)  =  ±^sin  u. 

It  has  been  proven  in  (i)  that  Ish^ul  >  l?l  Isin  ul  for  ^  €  (-«,») 

and  ^  *  0,  hence,  on  the  real  line,  (4.13)  is  satisfied  only  at  the  origin. 

Lemma  A. 2  hzis  shovnri  that  C  =  0  is  the  eigenvalue  with  multiplicity  of  2, 

—  T 

the  corresponding  cheu'acteristic  subspace  is  spainned  by  =  (cos  0,  -sin  0) 

_  T 

eind  =  (sin  0,  cos  0)  .  □ 

If  h  #  0  is  selected  in  (-iCj,»Cj),  there  is  no  pole  of  5J(t))  for  the 
Dirichlet,  Neumann  and  mixed  problems  on  the  line  Im  t)  =  h.  As  a  conse¬ 
quence  of  Theorem  4.  1  zind  Lemma  4.4  we  have  Theorem  4.2  and  Theorem  4.3. 


Theorem  4. 2.  If  f  €  ®  =  0,1,  0  <  h  <  then  equation 

(4.4)  with  boundary  conditions  (4.5)  has  a  unique  solution  u  e  Hp(D)  and 


Ia|  S  2 
(4. 18) 


id“Si 


I^(D) 


5  C  Ilf II 


H^(D) 


Z 

^=0,  1 


IC^I 


l-t 


(D)‘ 


||G°1|  _  (resp.  IlG^I  ,  ) 

H^(D)  H^(D) 


is  absent  In  (4.18)  for  the  Neumann  problem 


(resp.  the  Dirichlet  problem). 


Theorem  4.3. 
aind  then  the 
solution  u 

(4. 19) 


If  f  e  Lp(Q),  G"  e  i  =  0, 1,  0  <  8  <  1,  8  >  1  - 

equation  (4.1)  with  boundary  condition  (4.3)  has  the  unique 

in  w5(Q).  and 
P 
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where  llG^ll  ~  (resp.  HG^ll  )  drops  out  in  (4.19)  for  the  Neumann 

W^{Q)  W^(Q) 

problem  (resp.  the  Dirichlet  problem). 

Remark  4. 1.  The  proof  of  Theorem  4.2  amd  4.3  is  similau'  to  those  for  the 

Poisson  equation  in  [2].  The  shift  theorem  in  the  space  (Q)  for  general 

elliptic  systems  was  given  in  [21],  where  the  Agramovich-Vishik’ s  conditions 

are  used,  and  the  shift  theorem  in  (Q)  for  the  plane  elaisticity  problem 

Pt  P 

was  given  in  [26]  without  verifying  the  conditions. 


Remark  4.2.  The  shift  theorem  can  be  easily  generalized  to  amy  k  >  2,  namely 


-JlQ) 


s  c  [if II 


*  I 

£=0.1  0 


provided  f  e  ^(Q)  ,  ^(Q)  and  there  is  no  pole  of  3?(t))  on  the 

line  Im  Ti  =  k  +  1  -  6.  But  in  practical  problems  the  shift  theorem  is  less 

-  -£ 

applicable  for  k  >  2  because  of  very  strong  conditions  on  f  amd  G  in 

-  -I 

the  neighborhood  of  the  origin,  for  instamce  when  f  amd  G  aire  analytic 

auid  vanish  rapidly  at  to,  but  do  not  belong  to  W^(Q),  k  i  2.  Hence  the  shift 

P 

theorem  for  k  >  2  is  not  directly  applicable  to  these  problems  with  analytic 

or  piecewise  auialytic  data  (see  [22]).  For  this  reason  we  addressed  the  shift 

theorem  in  the  space  W^(Q)  only  for  k  =  2,  and  will  address  it  in  the  space 

p 

^(Q)  for  k  >  2  in  Section  5  . 

P 

Corol  lau*v  4.  1.  If  G^  vanishes  at  the  origin  and  G^€  H?  ^(Q)  then 

P 

(4. 19)  cam  be  rewritten  ais 


(4. 19' ) 


Hull  „  SC 

l^(Q) 


£»0. 1  "3 


J2-1  2-1  ] 

HX  {Q)J 


Proof.  By  Lemma  A,  2  of  [2]  G^  c  £  *  0,1,  and 

P 


-1 

1^11 


,  <  C  II  Gil 

VC  (Q)  IT  ^ 


3 


(Q) 


Then  (4.19')  follows  eaisily  from  (4.19).  □ 

Corollary  4. 2.  Let  u  =  A  denote  the  displacement  in  Cairtesian 

coordinates  and  u  €  V*5(Q),  0  <  0  <  1  be  the  solution  of  (4.1)  and  (4.3). 

P 

Then  u  €  vr!(Q).  aind  there  are  constants  c,  and  c„  >  0  independent  of  u 
p  Id 

auid  u  such  that 


(4.20) 


cJIull  _  <  Hull  -  <  C-  Hull  _ 

V^(Q)  V^(Q)  ^  V^(Q) 


Proof.  Since 


u,  =  u  cos  8  -  u„  sin  0 

1  r  8 

u-  =  u  sin  0  +  u„  cos  0 

2  r  0 


obviously  Hr^  ^ufl  _  *  Hr^  ^uH  ,  <  Hull  -  .  Further  note  that 

L‘'(Q)  L'^IQ)  VC(Q) 

p 


=.2  2  -2 

1  ^^1,1  ^ 

r^  ae^  r^  la0^ 


a^u. 


cos  0  - 


00 


1  2  r^'^r  ^^0  1 

®J  •  ^  ®  ^  ®J 


-X  (u  cos  0  +  u  sin  0) 

6  I*  U 


Then  we  have  Hr' 


0  £  k  i  2, 


3-2 


—  2 
I _ II  ^ 


^11  p  5  C  ||u||~  .  Similarly  we  can  prove  that  for 

ae^  L^(Q)  vC(Q) 

p 


l^p-2+kjjky  |2  ^  ^  II  ^11 2 

L‘(Q)  U^(Q) 


which  yields 


Hull  s  c-iiuii  „  . 

VIJ(Q)  W^(Q) 


In  the  exaictly  same  way  we  can  prove  the  other  half  of  (4.20). 
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For  the  regulzu'ity  of  the  wesik  solution  of  (3.1)  we  are  interested  in  the 

I 

special  auxiliary  problem  for  which  f  and  G  vanish  for  r  >  1. 

Theorem  4.4.  Let  u  €  W?(Q)  be  the  solution  of  (4.1)  and  (4.3)  with 

p 

f  €  L„(Q),  €  H?~^*^“^(Q).  £  =  0,1  0  <  ^  <  1.  3  >  1  -  k,,  and  let  u  = 

p  p  i 

_  -i 

be  the  displacement  in  Cartesian  Coordinates.  If  f  and  G  ,  £=0,1 
vanish  for  r  >  1,  auid  G*^  vanishes  for  r  =  0,  then  for  the  problem  with  the 
Dirichlet  boundairy  condition  (4.3a)  amd  mixed  boundau'y  condition  (4.3c) 

(4.21a)  +  l|r‘^5|l^  <  «  . 

L^(Q)  L^(Q) 

(4.21b)  I12)^u||^  +  llr'^ull^-  <  od  . 

L^(Q)  L^(Q) 

amd  for  the  problem  with  the  Neumann  boundamy  condition  (4.3b) 

(4.22)  <  «  . 

L^(Q) 


Proof.  We  first  prove  (4.21)  for  the  Dirichlet  problem.  Select  0'  such 

that  1  <  <  1  +  Kj.  Then  0'  >  1  >  0  amd  -Xj  <l-^'  <O<1-0<Kj. 

Due  to  Lemma  4.4  ?l(i?)  has  no  pxjle  between  the  line  Im  t)  =  1  -  ^  and  the 

line  Im  7)  *  1  -0'.  Because  f  and  G^,  £  =  0, 1  vanish  for  r  >  1  and 

G-  =  0  at  the  origin,  f  e  (Q)  amd  G^  e  wf,  fo),  £  =  0,1.  By  the  amgument 

u  p  p 

of  Theorem  4.3,  u  e  W?, (Q).  Thus  we  have 


(4.23) 


0; 


rdrde 


fXi 


r2(B'-2)|5|2l 


rdrde 


1^(Q) 


On  other  hand  we  have 
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(4.24) 


rdrd0 


S  J  I  »  r"<3-^'l5|2]rdrde  S  Il5l|2 


1^,(0) 


which  together  with  (4.23)  yields  (4.21a).  Noting  that 


,2  ,ml-i2  ,  -l-,2  2 

12)  u|  *  12)  u|  +  |r  u|  +  -; 


,  -1  ,2  ,  -l-,2 

Ir  u|  =  |r  u| 


“r  -  “e)  ^ 


we  have  (4.21b)  inunediately. 

The  proof  for  the  mixed  problem  is  similar  and  will  not  again  be 
elaborated  here.  Next  we  prove  (4.22)  for  the  Neumajin  problem.  Let  3'  be 
selected  as  above,  then  by  Lemma  4.4  the  origin  is  the  only  pole  of  XCtj) 
between  the  line  Im  tj  *  1  -  0  and  the  line  Im  tj  *  1  -0'.  Let  u  = 


u  (to  0)  and 
r 


.♦1(1-0')  . 

U*(T.0)  =  I  K(7,)(f.gh  e^’’’^d7j 

V§i  ''-c»+i{l-0' ) 


Arguing  as  above  u*  is  a  solution  of  the  Dirichlet  problem  (4.1)  and  (4.3b) 

in  W?, (Q).  By  Theorem  2  of  [26] 

P 

S  -  5*  ♦  ajSj  ♦  a2?2 

where  auvd  are  defined  as  in  Lemma  4.4.  Let  u  =»  A  and 

u*  *  A  ^u*,  then 

T 

U  *  U*  +  (Cj.Cg) 
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Note  that 


(4.25) 


■^0  "^0 


rdrd0 


=  J  J  ^|2)^u|^  +  Ir^ul^jrdrdO 


^  f  r 
•'0  "^0  ■ 


-,2  20-2)  ,-,2l  ^ 

u|  +  r  u  rdrdG 


^  l|u||\ 

W^(Q) 


Cl)  00  (l)  00 

(4.26)  f  f  ID^ul^  rdrde  =  f  f  \D^u*\^  rdrd0 
JnJi  JnJf 


S-l"  .  r"'«'-2>|S.|2 


I  u* I  rdrd0 


S  llu*||^ 


Combining  (4.25)  and  (4.26)  we  get  (4.22).  □ 

Since  |2)^u|^  =  |D^u|^,  |2)^u|^  =  iD^ul^  we  have  the  following  corollary. 

Corol  lairy  4. 3.  Let  u  aind  u  be  the  saune  ais  in  Theorem  4.4.  Then  for  the 
problem  with  Dirichlet  condition  (4.3a)  amd  mixed  condition  (4.3c) 


(4.27a) 


ID^ull^-  +  llr'^ull^  <  a> 
L‘^(Q)  L^(Q) 


(4.27b) 


ID^ull^  >  llr’^ulll^  <  . 

L‘^(Q)  L‘^(Q) 


and  for  the  problems  with  Neumann  condition  (4.3b) 


(4.28) 


IID^uM  <  «. 

L‘^(Q) 


Remark  4.3.  (4.21)  and  (4.27)  may  not  hold  for  the  problem  with  the  Neumann 

condition,  (4.22)  and  (4.28)  may  not  hold  for  u. 
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§5.  REGULARITY  OF  SOLUTION  OF  THE  ELASTICITY  PROBLEM 

We  shall  first  discuss  the  regularity  of  the  weak  solution  of  the 

problem  (3.1)  over  the  polygonal  domain  n  with  f  e  LAil)  and 

P 

t  JZ-t  2-1 

G  €  HI  ’  (Q),  namely  the  relationship  between  the  wealc  solution  and  the 

P 

solution  of  the  auxiliary  problem  over  an  infinite  sector  Q  with  f  and 

I 

G  having  a  bounded  support.  Then  we  shall  derive  the  regularity  of  the 

solution  of  the  problem  (3.1)  in  the  countable  weighted  Sobolev  space. 

Assume  that  A^  is  located  at  the  origin  emd  lies  in  positive 

x,-axis.  Let  (r,0)  denote  polar  coordinates,  and  let  5  6  (0,1)  and  S_  = 
i  o 

{(r,0)lO  <  r  <  5,  0  >  0  <  <j>  c  Q.  Let  ^-(r)  be  a  cut-off  function  in 

o 

00  1  5 

C  (R  )  such  that  s  1  for  0  <  r  <  =  and  s  0  for  r  >  5.  If  T,  ur„ 

o  do  1  M 

T  0 

c  r^,  let  V  =  (Vj,V2)  =  0^(r)u,  otherwise  v  =  ^g(r)(u  -  G  (A^)  where  u  = 

T  1 

(u,,u-)  €  H  (n)  is  the  weak  solution  of  (3.1).  v  is  extended  to  Q  =  S  = 

id  00 

{(r,0)|O  <  r  <  C9,  0  <  0  <  w>  by  zero  extension  outside  S_. 

o 


Theorem  5. 1.  Let  u  be  the  weak  solution  of  the  problem  (3.1)  with 

f  €  L-(n)  and  G^  €  H?  ^(G),  0  <  ^  >  1,  p  >  1-k,  ,  emd  let  v  =  if 

P  p  1  o 

Tj  vj  c  and  v  =  ^^(u  ~  G^(Aj)  otherwise.  Then 

(i)  If  Fj  c  r°  or  c  r°,  or  Tj  u  c  F°,  then  v  =  W^(Q). 


and 


(ii)  If  F,  w  F„  c  F^,  then  (v  -  V  c,e,)€U?(Q),  where  e ,  = 
1  M  2  i  P  1 

T  T 

(1,0)  ,  *2  “  .  Cjt  i  *  1.2  are  some  constants,  and 


(5.1b)  Mv  -  Y.  =i«ill.^_,  ^  C  ^||f||j_ 


i»l,2 


W^(Q) 


(S  )  ^  “  1  1  ^  1 

»  «  (S3)  h1(S3\S3/2 


J 


24 


Proof.  (i)  We  prove  (5.1a)  for  the  case  F  u  c  F  (Dirichlet  boundary 

1  N 

type),  the  proof  for  the  mixed  boundary  type  is  similar  with  what  follows. 


At  first  we  may  sissume  that  G  =  0  at  then  v  satisfies 


(5.2) 


(J-Cv)  =  +  Lj(u,  ^g)  =  f  in  Q 

0=0,6)  ~  ^3^  1 0=0,  6)  “  ^  1  0=0,6) 


where  L^Cu.^^)  is  a  sum  of  terms  D^uD**  with  0  <  |a|  <  1,  1  <  |a' I  ^  2 
Obviously  f,  and  v  vanish  for  r  >  3  and 


(5.3) 


^"l  (Q)  -  ^  (S  )  ^  1  1 


For  w  e  H^(Q)  =  {w|  ||D^w(|  >  “Iq  r>  =0)  we  have 

°  L^(Q) 


B(v 


,w)  =  f 


f  •  w  rdrdir  ,  Vw  €  H^CQ) 


where  B  is  the  bilinear  form  defined  in  (3.4). 

On  the  other  hand,  by  Theorem  4.3  and  Corolleu'y  4.1~4.2  there  exists  a 
unique  solution  z  =  of  (5.2)  in  w|(Q),  and 


(5.4) 


-^(Q) 


h!(Q)' 


3 


S  C  I  Ilf  II  _  .  +  IIG'  II  +  ||u„ 


.1 


)■ 


and  owing  to  Corollary  4.3  ||D  zH  _  <  oo  . 

L^(Q) 

Let  Hq(Q)  =  {w|w  €  Hq(Q)  with  bounded  support  in  Q} .  Thus  we  have 


b(2,w)  =  r 

Jn 


f  •  wrdrd0  ,  Vw  €  Hq(Q); 


hence 
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B(z-v, w)  =  0 


Vw  €  Hq(Q). 


Since  H^(Q)  is  dense  in  Hq(Q),  B(z-v, z-v)  =  0.  Therefore  (z-v)  can  only 

be  a  rigid  body  motion.  Due  to  the  Dirichlet  boundary  condition 

V  =  z  €  wJCQ). 

P 

If  *  0  at  the  origin,  let  u  =  u  -  amd  v  =  Then  we 


have 


L(v)  =  +  Lj(u. =  f  in  Q 


''le=o.»= 


=  G 


6=0,  (j 


0=0,  Cl) 


Obviously  f  €  L^CQ)  ,  f  and  G  vanishes  in  Q  for  r  >  5.  Applying  the 
P 

results  above  to  v  we  have  that  v  =  ^-(u  -  G*^(A, ))  e  W?(Q)  amd  (5.1a) 

o  ip 

holds. 

(li)  If  r,  u  Fw  c  r\  v  satisfies 
i  n 


(5.5) 


L(v)  =  +  L^(u,0g)  =  f  in  Q 

T(v)|  =^g^+^(a,^)|  =gl=G^| 

0=0 ,  u  0=0 ,  w  0=0 ,  u) 


where  Lj^(u, ^^),  f  ame  the  same  ais  in  (i),  and  ij(u, ^g)  consists  of  terms 

uD%.,  |a|  =  1.  Obviously  G^€  W^CQ),  amd 
o  p 

(5.6)  IIG^I  -  5  C  fllG^I  -  ,  +  Hull  ,  1  , 

V^(Q)  H^''(S^)  h'(S3\S3/2)J 

and  f,G^  amd  v  vanish  for  r  >  5.  For  any  W  e  H^((3)  =  (WIUD^WH  _  <  «> 

L^(Q) 

B(v,w}  =  I*  f'wdx  ♦  f  g^'Wds. 

•'q  '^aq 

On  the  other  hand,  by  Theorem  4.3  and  Corollamy  4. 1~4.2  there  is  a  unique 
solution  z  ■  (Zj,Z2)^  of  (5.5)  in  amd 
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(5.7) 


Izll  _  <  C 


fllfllL  (0,  '  l£‘l  1  1 

^  wi{Q)J 

P 


by  (5.3)  and  (5.6) 


S  C 


(./II 


,  (S  )  *■  1  I  *  1 

S  «’  (Q)  « 


]• 


zmd  by  Theorem  4.4  HD  zH  _  <  «. 

L^(Q) 

Let  H^(Q)  =  (wlw  e  H^(Q)  with  bounded  support  in  Q> ,  then  z 


satisfies 


B(z 


fwdx  +  I  g^‘v#ds,  Vw  e  H^(Q). 


,  w)  =  r  fwdx  +  j  1 

••q  ''ao 


a, 


Therefore 

B(2  -  v.w)  =0.  Vw  €  H*(Q). 
Since  H^(Q)  is  dense  in  H^(Q)  we  have 


B(2  -  V,  2  -  v)  =  0, 

which  indicates  that  the  strain  energy  is  zero,  z  -  v  represents  only  a 
rigid  body  motion,  i.e., 

3 

'  -  ’  =  Z  'I'l 

i=l 

T  1 

where  =  (-y,x)  ,  and  Cj^,  i  =  1,2,3  au'e  constants.  Because  D^(z  -  v) 

2 

€  L  (Q),  c^  =  0.  Then  (5.1b)  follows  from  (5.7). 


Corollary  5.1. 
Proof.  If  Fj 
Rq  >  0.  Hence 

in  Q.  Since 


V  is  continuous  in  Q,  auid  (v  -  v(A, ))  e  V?(Q). 

1  p 

1  ^  ^  2 
u  Fj^  c  r  ,  by  Theorem  5.1.  (v  -  £  ^^e^)  €  (Sj^)  for  any 

2  0 

(v  -  £  c.e, )  €  C  (S_  )  (see  [2]),  therefore  v  is  continuous 

1.1  1  I  Ro 

2  2 
(v  -  £  c  e  )  e  H?(Q),  0  <  ^  <  1,  r^"^(v  -  £  c.e.)  e  L^(S_  ) 

1.1  1  1  P  1,1  1  i  Ro 


27 


Hence 


2 

which  implies  that  J]  c.e. 

i=l 


v(A, )  and  that 

X 


(v  -  v(Aj))  e  w|(Q). 


If  Fj  c  r°  or 

V  e  )  for  any 

P  Ho 

continuous  in  Q  eind 


r„  c  r°  by  Theorem  5.1  v  =  ^.(u  -  G°)  e  W?(Q). 

M  op 

R_  >  0,  thus  V  -  G^(A. )  €  ),  and  v  is 

u  1  Kq 

v(A^)  =  G*^(Aj),  and  (v  -  v((Aj))  €  W^(Q). 


T  —  T 

Corollairv  5.2.  Let  u  =  (u,  ,u_)  and  u  =  (u  ,u  )  be  the  weak  solution  of 

X  ^  Vo 

(4.1)  in  Cau'tesieui  and  polar  coordinates,  respectively.  Then  (u  -  u(Aj^)) 

1 


(5.8a)| 


u-u(Aj^)| 


H^(S 

S/2 


^  S  cj^iifii 


•Z 

1=0 


S,) 


]• 


1 


(5.8b)||5-5(Aj)| 


W^(S 

3/2 


^  ^  c[^||f|| 


■I 

t=0 


IG^I 


U 


,]■ 


Proof.  Note  that  v  *  u  and  v  =  u  in  then  (5.8a)  and  (5.8b)  follow 

easily  from  Theorem  5.1  and  Corollary  5.1. 

Remark  5. 1.  The  regularity  of  the  weak  solution  of  the  elasticity  problem  in 
polynomial  domain  was  addressed  in  [261 .  It  was  concluded  in  [26]  that  v  = 

V  ‘  '^(0)  with  0=1+  e,  c>0  arbitrary  provided  f  € 
t  2.-1  2-1 

G  €  w  ’  (Q),  and  K(ti)  han  no  pole  on  the  line  Im  t)  =  -e.  Actually  the 

condition  0  >  1  is  not  necessany.  Theorem  5.1  indicates  that  v  e  H?’^(Q) 

0 

or  (v  -  v(Aj))  €  w|(Q)  if  0  >  1  -  Kj.  (Cj  is  the  smallest  positive 
imaginary  part  of  eigenvalues  of  the  operator  3{(tj),  which  depends  on  the 
geometry  of  the  domain,  the  type  of  boundary  conditions  and  the  material 
properties.  The  condition  0  >  1  -  ic^  precisely  reflects  the  nature  of  the 
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singulau'ity  of  v.  If  k.  <  1,  then  v  e  H?’^(Q)  with  3  €  (1  -  k,,l),  or  if 

1  p  1 

kj  >  1  then  v  €  and  V  can  be  even  smoother.  We  will  elaborate  on 

this  later.  Nevertheless  there  always  exists  some  3  ^  (0,1)  such  that  v  e 

H?’^(Q)  and  (v  -  v(A, ))  e  U?(Q).  The  improvement  above  is  substantial. 

3  13 

Lemma  5.  1.  For  3€(0,1),  k>^2:0,  0  ^  i  <  2  and  u(Aj^)  =  0  there  exists 


Cl,  C2  >  0  such  that 


(S.9) 


cJIull  . 

K3  ^5^ 


-  ic  £  -  k  £ 


Moreover,  if  for  some  constant  C,  d  >  1  zuid  t  <  |a|  ^  k 


(5.  10a) 


|r“^"^‘"3D“u||  ,  <  Cd^‘^(k-£)! 

l2(S^) 


then  for  |a|  =  k  and  C  ^  MC,  d  <  Nd, 


(5, 10b) 


^  Cd^*^(k--£)! 

^  'S3) 


Vice  versa,  if  (5. 10b)  holds  for  £  ^  |a|  £  k,  then  (5.  10a)  stands  for 
|a|  =  k  with  C  5  MC  and  d  S  Nd. 


Proof.  Note  that 


k  k 

a  u  a  u^ 

=  cos  0  — n - sin0  — ^ 

ar*'  dr 


k  k 

a  u  a  u^ 

sin0  — ^  +  cos  0  — j— 

ar*  ar 


Then 


(5. 11) 


k-^+3  5  u,, 


k-<+3  d\, 


k"  2 

ar*^ 


k"  ? 

ar*^  L‘^(Sg) 


s  11  ul 


'  k  I 


which  is  the  second  inequality  of  (5.9)  for  *  I  “I  ®  for*  “2  ~ 


Ia|  =  k 


we  have 


Therefore 
(5. 12) 


m=0 


■  ,k-m 
d  COS0 

-m 

3  u 

r 

d^  ®sin0 

a®Ug. 

.-k-m 
■  d0 

30® 

de*'-" 

30®  ■ 

■  ,k-m  . 

3®u 

jk-m 

Q® 

a  u„- 

d  sin0 

d  COS0 

0 

■  de‘‘‘" 

as" 

d0 

30®  ■ 

30^  L^(sj  ^  W  ae® 


since  0  <  3  <  1  and  uCA^)  =0  ,  due  to  Lenuna  A. 2  of  [2]  we  have 


Ir^‘^u|| 


Ir^'^ull 


2  -  1  1 
L^CS^)  ° 


2  ^  ^0"“*'  1  1 

L^CS^)  Kg-'is,) 


s-2  au,, 


s  c^lluH 


L  (S3) 


Hence  for  0  ^  I  ^  2,  k^£,  we  have 

^  a  ^k _ 


30%  2,^  ,  -  '-0'"*"  2,2 


(S3) 


(5. 13) 


30^  L^(S,)  ^  W  30®  L^{S,) 


5  C  Hull  .  , 

(S5) 

Actually  we  can  make  similar  arguments  for  each  term  of  type  r*^  0  < 

5  |a|  *  k.  Hence  we  have  the  second  Inequality  of  (5.9).  The  first 
ineqxiality  of  (5.9)  can  be  proved  in  exactly  the  same  way. 

Next  we  shall  prove  (5.10b)  if  (5.10a}  holds.  (5.10a)  and  (5.11)  ~ 
(5.13)  lead  to 


(5. 14a) 


^  Cd^-^k-«! 

ar*  L‘^(S3) 


30 


(5. 14b) 


m=£ 

C^l+ij  (l+d)^"^(lc-£)! 


which  au'e  (5.10b)  for  =  0  amd  a.^  ~  |a|  =  k,  respectively.  Similar 
estimates  can  be  derived  for  general  terms  r*^  0  <  ^  la|  =  k. 

Hence  (5. 10b)  holds  with  d  ^  Md  and  C  5  NC  where  M  =  2  and  N  =  C_C 


('  *  r- 


In  exactly  the  same  way  we  can  prove  that  (5. 10a)  holds  for  |a|  = 


k  with  d  i  Md  and  C  ^  NC  if  (S. 10b)  holds  for  £  =  |a|  <  k. 


Theorem 


5.2.  Let  f  e  H^’°(fl).  g^  *  €  H^'^^“^n),  £  =  0, 1  p  = 


0^.02' •••  0  <  <  1,  >  1  -  k|  (Kj  is  defined  in  (4.15)  with  respect 

to  the  vertex  Aj^),  1  ^  i  2  M,  and  IF^I  *  0.  Then  the  problem  of  (3.1)  has 


unique  solution  u  €  H^'^^’^(n),  and 

p 


(5. 15) 


H^"'^'^(n) 

P 


5  cjllfl 


„k.o,n)  L.  ”  V 

**3  £=0,1  "p 


Furthermore,  if  f  e  B°(n,C_,d.),  G^  €  B^"^n,C  ,,d  ,),  £  =  0,  1  then 

P  f  f  P  j,£  g£ 


u  €  B^(Q, C  ,d  )  with  C  and  d  satisfying 

p  U  U  U  U  ^  'ey 


Cu  i  (Cj  >  Cgo  ^  C^i). 


(5. 16) 


du  i  max(dj.,djjp.clj,j) 


where  and  aure  some  constants  independent  of  f  and  G^ 


>.2.  If  ir^l  »  0,  the  theorem  holds  provided  f  and  satisfy 


usual  condition: 


r  f  dx  +  r  gjds  =  0,  1  =  1.2 

Jjj  1  Jj.  1 

JjjVl  ■  Vl’*"  * 

Then  the  solution  uniquely  exists  (up  to  a  rigid  body  motion)  in 

p 

and  B?(n).  respectively. 

p 

Proof.  For  the  sake  of  simplicity,  let  Q  be  a  straight  line  polygon  shown 
in  Fig.  2.1.  By  Theorem  3.1.  the  problem  (3.1)  has  the  unique  solution  u  = 
(Uj.U2)^  e  H^(n). 

Let  S,  _  =  ((r,.0  )|0  <  r,  <3.,  0  <  0.  <  w.}  c  n.  shown  in  Fig.  5.1. 

l.oi  li  11  ii 

where  rj^  and  0^^  are  the  polar  coordinates  with  respect  to  the  vertex  Aj^ 
and  edge  F^.  Assume  that  0  <  <  1  such  that 

^i.23  ^  ^J.28  *  ^  i  J*  *  1*2 . M. 


Vi/ 

\G\  1 


Fig.  5. 1  Neighborhood  of  Vertex 


Let  Q,  »  n\  U  S 


L  »  n\  U  S,  _  .  By  the  au'gument  of  difference  quotient  it  is  easy 
i*l  ^’*^1 


to  show  that 


(5. 17) 


'“V‘==(n  ,  '  ‘'("'VcQ  ,  *  ^  )] 

"  ^“8/2^  ^”s/4^  «»0,1  "  ^^8/4^ 


32 


5  C^l 


Jc  0  ^ 

(Q) 

p 


X  ’*^+3/2-^.3/2-£  » 
£=0,1  3  (r*^)‘ 


If  f  €  Bp(n),  6  ^(n),  l  =  0,1,  by  Morrey’s  argument  in  '^231  we  can 

prove  that  u  is  analytic  in  Q  and  on  F  except  the  vertices,  suid  there 
exist  some  constants  Cq  and  d^  satisfying 


(5. 18) 


S  ^  «0  *  =  0  *  =  1> 

G  G 


do  ^  Nq  max{d^,d  d  ) 
G  G 


such  that  for  |al=ki2,  l^i^M 


(5. 19) 


la 

■-  'Va' 


By  the  arguments  of  Theorem  2. 1  and  Lemma  5. 1  we  have 


(5.20) 


"•1“"^"'  2,,  ,,  ,  ^  Vo"""'  - 


where  Cq  and  dQ  may  be  different  from  those  in  (5.19),  but  we  use  the  saune 
notations  for  simplicity. 

By  Theorem  2.1  auid  Lemma  5.1  f,f  6)Co’^(n)  (resp.  »2(n)),  g^.g^  e 

P  P 

„k+3/2-£,  3/2-£,_£..  ,  „3/2-£,_£,,  t-r-t  t  t.  *.  t-u  i.  • 

)€„  (r  )  (resp.  ®  (r  )),  and  it  is  sufficient  to  prove  that  in 

P 

each  sector  S,  ,  ,,,,  1  5  i  5  M 
1 , 


^8  ^^i,3^/2^  ^8  ^^3^^  £=0,1 

+  Hull  1,^1  ] 


k+2,£,  2-£ 


where  depends  on  k,  aind 


(5.22) 


^  ^  L  fS  )  111 

respectively,  where  |a|  =  k  +  2,  k  i  0,  L^,  Pj^  amd  Dj^  are  sufficiently 
large  constants,  but  Independent  of  k.  ia^  ~  2)  =  -  2  if  ^  2  and 

-2  =  0  if  <*2  <  2. 

By  Corollary  5.2  we  have 


(5.23)  ||u  -  u(A  )|| 

W\a.y2^ 


(S  )  "  I]  L2-l,2-l 

^  1^0,1  f^i.d. 


«  (Si.5.^Si.5^/2^- 


We  may  aissume  that  i  =  1,  is  at  the  origin,  lies  on  the  positive 

Xj-direction,  and  assume  without  losing  generality  that  “(Aj)  =  0.  To 


simplify  the  notation  we  shall  write  S,  _  =  S-,  B, 

1,  o  j  o  i 

three  cases  to  be  considered: 


=  B,  etc.  There  are 


(i) 

r„  c  r°. 

(ii) 

^1* 

r„.  c  r‘. 

(iii) 

c  r°.  c 

In  case  ( i ) 

we 

assume  G  = 

-k  kd^u 


V  satisfies 


(5.24) 


k”2,  By,  ,  _ 

L(v  )  =  r  (r  f)  in  S, 


^  ®  ' 
0*0,  u 


Then  applying  (5.23)  to  the  problem  (5.24)  we  obtain 


34 


(5.25)  Hub'll  < 


5 


=  f,,  k-2a^  ,  2»,,,  ,,  k  a^u  „  , 

r  3^1=  ’"S'=S>  ar''  K  ‘VS3/2>. 


which  Implies  that  for  |a|  =k+2,  k>0,  0  i  i  2 


(5.26)  2)“u||f  .  1  ^  cfr 

5/2^  1“ 


k-2  a^(r^f)„  , 

ar*^  Sr^'  ^  ^^5'^^5/2^' 


-<c[ 


llfll  ^  0  >  Hull 

<  '=5)  K  'V®S/2>' 


For  ttg  >  2.  for  instajice,  a2  =  k  +  2,  k  >  0  we  have  from  the  equation  (4.1) 


(5.27)  -p 


.k+2  _k.  -k+2  .k+1  _k 

au  ^3  t  /-^au  au  au 

r2r,.,^vf2  r  r  r 

-kV2  "  —r  *  r  ““2~k  - k  — r 

as  ao^  t  an  ao  arae*^  ae^ 


a  u  .. 


ak+l 

3  u„ 


.k+2 

a  u 


+  (X  +  3p)  — j— y-  +  (A  +  p)  r  - 

ae  arae 


e _ 

k+l  ’ 


-k+2  .k_  .k+2  ^k+1  .k 

,=  ,0.  ,  >  8  2  0  r2®  "9  ^  “9  ®  “9l 

(5.28)  (X  *  2|i)  — jjjj-  .  r  —  *  |l  k  ;  k  ^  '' - k - iT 

ae  ae  '•  ar  ae  arae  ae  J 


-k+l  _k+2 

a  u  a  u 

+  +  3p)  ■  /  +  (A  +  m)  r  - ^ 

ae  arae 


(5.27)  amd  (5.28)  lead  to 


,k+2- 


(5.29) 


gk+2- 

)  '*'  II  5  k  11 L  (S  ) 
>  ar  ae  6^ 


arae*' 


-2  a^'u 


ae*' 
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o  a*'"*’!- 

II  ^  U,,  „  “ 


'  1 
araa"** 


with  C  -  maix 


'  •  a- 


Suppose  (5.21)  holds  with  u  =  0  for  OS  |a|  S  k  +  2  and 
0  S  *2  S  k  +  1.  Then  we  have 


,k+2- 


(5.30) 


ae''*='S'®a/2'  ' 


11  V  n  k+i 

K@  'Sjl  K*  (Ss''S3/2)- 


which  is  (5.26)  for  ag  =  |al  =  k  +  2.  Actually  we  can  prove  (5.27)  for  each 
term  of  type  r*^  ^2)*u  in  exactly  the  saune  way.  If  *  0,  let  v  =  u  - 
then  applying  the  results  above  to  v  we  get  (5.21)  with  absence  of  G^. 

Next  we  shall  prove  (5.22)  for  case  (i).  Assume  again  that  G^  =  0,  let 
f  €  B°(n.  Cp  djf)  then 

fB  II  k-2a*'(r^f)  ,,  ^  ,,  ka^f,,  k-ia*'"^?  ,, 

(S.31)  llr  llr  Sklr  ^ 


♦  k(k  -  1)  llr 


k-2  , 

ar‘^-2"s^^a^ 


S  C,[d*  *  2d^-'  *  d^-2]k! 


S  3CjdJk! 


Also  due  to  (5.20)  we  have 


(5.32) 


kaV 


ak+1- 

d  u„ 


ar*' ar*'  '  L^(S,\S,,^) 


flk- 

♦  n~\\  2 


"  “F  2 

ar  ae  l^(s^\S3/2) 
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S  CQl^d^'^dc-l)!  +  -  2)!  +  cIq  ^(k  -1)! 

=!  -  »>' 

Substituting  (5. 31 )~(5. 32)  into  (5.26)  we  get  (5.22)  for  |a|  =  k  +  2.  k  >  0, 

0  i  i  2  with  L  =  4C(Cq  +  C|:)  and  D  =  max  (d^,d_Q). 

G 

Now  we  shall  prove  (5.27)  for  a2=k+2,  k>0.  Suppose  (5.22)  holds 
for  OS  |a|  S  k  +  2  and  0  <  a2  -  ^  (5-29)  we  have 

(5.33)  1  ^  clc^Jfk!  +  LoV'^k!  +  LD*'‘^P^‘^(k  -  1)! 


■f  LD^'^P^  ^(k  -  2)!  >  LD^*^P^"^(k  -  1 ) !  ■*•  LD^P^ 


S  LoVk!  ^ 

L  lUKi 


_C _  . 

P^Dk(k-l) 


k  k 

S  LD  P  k! 

where  P  and  D  are  selected  large  enough,  e.  g.  P  =  4C  auid  d^/DP  <  1. 
Similau'ly,  we  can  prove  (5.22)  for  each  term  of  type  r*^  2  S  a2  S  |a|  = 

k  +  2.  For  G°  *  0.  let  w  =  u  -  G°  and  f *  =  f  -  L(G°)  €B°(fl.  Cj:,,  d^,) 
with  Cjf,  =  ■*■  (5^  *  lit*)  CgO  and  dj*  =  max ( d , dgO ) .  Then  applying  the 

results  above  to  W  auid  f  we  get  (5.22)  in  general  (G^  is  absent),  with  L 

a  Mq(Cq  *  *  C_q)  and  D  *  max(d|;,d_Q,  d^). 

^  k-  ® 

In  case  (ii):  ^  =  r*'  ^  ,  k  >  2  satisfies 

ar* 
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(5.34) 


k-2  a^(r^f) 

L(v  )  =  r  — - 

ar 


in  S, 


'“■re-  “e* 


k-1  a^lro') 


l0=O,u 


ar  10=O.U 


Applying  (5.23)  to  (5.34)  we  have 


(5.35)  11^ 


^  r„_k-2  a*'(r^f)„  .  „_k-l  a^(rG^ )  „ 

I  5  5  ^  c  l||r  - j- — IL  .  +  l|r  - 1 

I  ar"  ar  K^’  (s^) 


-k 

+  V 


} 


which  implies  that  for  |a|  =  k  +  2,  k  >  0,  0  ^  ^  2 


/r.  ,._at~2_a-„  ^  p;  f„_k-2  a^(r^f )  „ 

(S.36)  Hr  I  -  (I'-  'l^(S5 


^f),,  ^  ,,  k-1  a^(rG^)„ 

1  a^“ 

.  „_ka  u„ 


(llfl  t  n  *■  k*l  0  *  k*l 


'3  3 


V®a/2' 


^^5'^^3/2^‘ 


For  2  <  «2  ^  1^1  =  k  +  2,  arguing  as  in  case  (i),  we  csm  conclude  (5.22)  with 
the  absence  of  G^  for  0  $  ttg  5  |a|  =  k  +  2. 

If  f  €B°(n,  C|:.  d^),  G^  (5.32)  holds,  and 


(5.37)  II 


k-1  a*'  ^(rG^)„  ,  „_ka^  V„  ^  „_k-iaV„ 

~kri — II  11  -  11^  ~k7i“"L  (S  )  11^  ~ir“L  (s  ) 

ar*'  ^  Kp  (S-)  w  5r 

P  O 


k-l=l 


*  k(k  -  1)  l■-'‘■^**^E=r'L.(S.)  * 


k-1  a*'‘^^G^ 


ar*'"'  ar*'a0 


,k-2  a*'G'  .. 

ar*'"'a0 
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- 


+  (2k+l)  cl5^~^(k  -  D!  +  cl^~^k! 


+  dT,  k!  +  d2 
G^  G 


S‘  H 


^  6  C  k! 

G^  G^ 

which  together  with  (5. 31 )~(S. 32)  and  (5. 35)~(5. 36)  yieid  (5.22)  for  |a|  =  k 

+  2,  k  >  0,  0  S  ^  2  with  L  =  6C(Cq  C_j)  and  D  =  max(d|r.  d_j, 

G  G 

d-).  As  arguing  as  in  caise  (i)  we  can  conclude  (5.22)  in  general  for  case 


(ii). 


The  proof  for  case  (iii)  is  similar  to  csises  (i)  and  (ii).  We  will  not 
repeat  it. 

Summarizing  the  analyses  in  each  sector  euid  interior  we  conclude  that 
u  €)t^'^^’^(Q)  and  (5.22)  holds  if  f  6  JC^’°(n)  and  G^  €  ^'^(Q) ,  ^  = 

0,1.  Furthermore,  if  f  eB°(Q)  and  G^eB^'fn),  ^  =  0,1,  then  ueBo(n) 

P  p  P 

with  C-  and  d-  satisfying 

C-  i  M(C«  +  C  +  C  +  C-) 

“  G^  G 

'*S  =  <*-1.  ’‘o> 

G  G 


which  together  with  (5.26a)  imply 


=5  "  "c“'f  *  V  * 

G  G 


“3  ^  *  V  "  V’- 

G  C 


By  Theorem  2.1  and  Lemma  5.2  u  e  B„(n,  C  ,  d_)  with  some  C  and  d 

p  u  f  u  u 

satisfying  (5.16). 
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Remark  5.3.  If  kJ  =  1  then  C2in  be  chosen  amy  positive  number.  If 

>  1,  the  solution  u  may  be  smoother,  for  instance,  l<(c|<2,  l<i^M 

and  f  e  (resp.  ^‘^(Q)  (resp.  with  B,  € 

P  p  p  pi 

(icj  -  1,  1)  and  k  >  1,  then  u  €  (resp.  B^Cn)).  In  general,  for 

1  p  p 

kJ  e  (n,  n+l],l  5  i  5  M.  u  e  (resp.  u  e  Bo'^"(n)  if  f  and  G^ 

X  p  p 

are  given  in  H^’^(n)  auid  n+2  (resp.  B”(n)  and  b”^^  ^(G)). 

P  P  PR 
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APPENDIX 


The  transcendental  equations  (4.11)  ~  (4.13)  are  given  in  Lemma  4.1 
without  proof.  We  shall  derive  these  equations  here.  Consider  the 
homogeneous  problem 

(A.l)  1/(tj)u  =  [UDq.t)).  B(Dg.  Tj)]u  =  0 

naunely 

d  u  du 

-n  — +  (2>1+A)(1+7I  )u  +  (iT)(/i+X)  +  (3p+A))  =  0 

(A.2)  .  in  I  =  (O.u) 

du  d  u 

-  ((34+A)  -  iT)(p+A))  — - (2fi+A)  — 5-  +  fid+T)  )u  =  0 

de 

with  one  of  the  following  boundary  conditions: 

A 

(A. 3a)  u|„  -  =0  (Dirichlet  condition) 

o~U,  0) 

(A.3b)  =  °  (Neumann  condition) 

(A. 3c)  *  0.  =  0  (Mixed  condition) 


where  <r  and  <r  are  given  by  (4.10). 

To  00 

Since  the  coefficients  of  the  differential  operator  L  and  boundary 

^  ^  bo 

operator  B  are  constaints,  the  solution  u  can  be  written  in  e  C  = 

bO  T 

e  ^  suid  C  satisfy  A(b)C  =  0  with 


A(b) 


-  b^ji  ♦  (1  +  T)^)(2p  +  X) 
(iT)(M  +  X)  -  (3m  +  X))b 


(i7j(M  +  X)  +  (3»x  ♦  X))b 
-  (2m  +  X)b^  +  m(1  +  Tj^) 


Then 

(A. 4)  det  A(b)  *  m(2m  +  X)((b^  -  (1  +  tj^))^  +  4b^]. 


Let  T)  *  iz,  by  solving  for  b  from  det(A(b))  *  0  we  have 

bj  =  i(z  +  1),  =  i(z  1),  bg  *  (i(z  -  1),  b^  =  -i(z  -  1)  . 
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For  z  *  0,  ±  1, 


A(bj)  =  (z  +  1) 


(A  +  3fi)  -  z(m  +  A) 
-i(z(ji  +  A)  +  {3fX  +  A)) 


i(-z(fi  +  A)  +  (3fi  +  A)) 
(A  +  3fi)  +  (fi  +  A)z 


IT  ^ 

the  corresponding  eigen  vector  C  *  (1,  1)  .  Similarly  we  have  u  = 

(1,  -i)^,  ®  (1,  iH)^,  =  (1,  -iH)^  corresponding  to  b^,  and  b^ 

with 


(A.5) 


H  -  z(p-*-A)  ->•  (3n+\) 
z(4+A)  -  (3p+A) 


Therefore,  the  solution  of  homogeneous  problem  has  the  following  form: 


(A.6) 


P  fcosCz  +  1)0  I  p  fsinCz  +  De]  .  p  f  cos(z  -  1)0 
“l^-sinCz  +  l)0j  “2{cos(z  +  DoJ  j[-Hsin(z  -  1)0 


sin(z  - 


4lHcos(z 


-  1)0'! 

-  l)0j  • 


If  z  »  ±1  (i.e.,  Tj  =  ±i),  bj  =  2i,  b^  =  -2i.  b^  =  b^  =  0.  ACb^)  =  A(b^) 
are  null  matrix.  The  corresponding  eigen  vectors  are 

Cj  =  (1.  i)^.  =  (1,  -i)^,  =  (1.  0)^,  =  (0.  1)^ 

and  the  solution  of  homogeneous  problem  hats  the  form 


^  „  f  cos20  1^0  rsin20  1  ^  n  fll  ^  n  fO') 

“  “  ®l[-Msin20j  ®2[Mcos2eJ  ®3[oJ  ®4[lJ 


with 


1  for  z  =  1 


2p+A  for  z  »  -1 


For  z  »  0,  bj  *  b^  =  i,  bg  =  b^  ®  -1,  the  corresponding  eigenvectors  are 


Cj  »  (1,  i)^.  Cg  *  (1.  -i)^. 
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and  the  solution  of  the  homogeneous  problem  has  the  form 


U.Q) 


COS01  Q  fsinfll  o  f  +  X)0cos0  +  Asin0'l 

'*  ~  ®l(^-sin0j  2[_cos0j  +  X)0sin0  -^icos0  J 


^  _  rOfx  +  X)0sin0  -  Xcos0l 
°4[(3fi  +  X)0cos0  -  |isln0j 


B  =  would  be  determined  according  to  the  types  of 

boundau'y  conditions. 


Lemma  A.  1.  Let  t)  *  iz  be  an  eigenvalue  of  the  operator  IKtj)  for  the 
Dirichlet  problem  (A. 2)  and  (A. 3a).  then  z  satisfies  the  equation 


(A.9) 


,  2 

sin  Z(i) 


,  2 

sin  (j. 


Proof.  For  the  Dirichlet  boundary  condition  u|_  _  =  0,  if  z  *  0,  ±1 

f  Cl) 

to  (A. 3a)  and  (A. 6),  B  satisfies  £B  »  0  with 


£ 


0 

cos(z-l)« 
-sin(z-l )u 


0 

1 

sinCz+l )u 

COS(z+l  )(i) 


1 

0 

cos(z+l  )w 
-Hsin(z+l)u 


0 

H 

sin(z+l)u 
Hcos(z+l )« 

J 


For  the  existence  of  non-trivi2i21 


B  it  is  sufficient  and  necessary  that 


(A. 10) 


det(£)  =  (1  +  H)^  sin^w  -  (1  -  H)^  sin^zw  =  0. 


Note  that  -  *  and 

4  l-2w 


(A.  11) 


f  1+-) 

1+H  z(*i+A)  __  it' 

1-H  *  3(m+X)  “  ^ 


(A. 10)  and  (A. 11)  yield  (A.9). 


z 

3-4i; 


due 
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If  z  =  ±1,  by  (A.7)  B  satisfies  ZB  =  0  with 


cos2u 

-Msin2« 


sin2<i> 

Mcos2(t> 


and  det(Z)  =  4M  sin  w.  Since  0  <  w  5  2ir,  for  w  =  ir,2ii  the  homogeneous 
problem  has  a  non-trivial  solution.  But  z  =  ±1  euid  u  =  n,2n  satisfy 
(A.  9).  Hence  z  =  ±1  are  the  zeroes  of  (A.  9)  for  u  =  7i,2n,  which  are 
included  in  the  equation  (A. 9). 

If  z  =  0,  by  (A. 5),  B  satisfies  ZB  =  0  with 


-sln« 


sinu  Ofi-t-AlcdcostiH-Asintt) 

cosw  -(3fi+X)(Jsinw-ficos« 


( 3n+\ ) ws i nw- Acosw 
( 3fi+ A )  ucosw-ps  i  nw 


2  2  2  2 

and  det(Z)  *  (3n+A)  «  -  (X+ji)  sin  «  >  0  which  implies  B  =  0.  Hence  zero 

is  not  an  eigenvalue. 


Lemma  A.  2.  Let  ij  =  iz  be  8ui  eigenvalue  of  the  operator  IKtj)  for  the 
Neumainn  problem  (A. 2)  ^Lnd  (A. 3b),  then  z  satisfies  the  equation 


(A. 12) 


2  2 
sin  zu  a  z  sinu. 


specially  t)  >  0  is  an  eigenvalue  with  multiplicity  of  2,  the  corresponding 

—  T  —  T 

eigenvectors  are  e^  »  (cos0.  -sinO)  aind  e^  ^  (sin0,  cos0)  . 


Proof.  For  Neumann  boundairy  condition  <r  I  *  <r  I 

(4. 10) 


*  0  we  have  by 
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(A. 13) 


*  -  I’^e)  ■  “ 


<r„-  =  (24+A)  +  (24+X+\z)u  =  0 

oo  aa  T 


For  z  *  0,  ±1  by  (A. 13)  and  (A. 6),  B  satisfies  ZB  =  0  with 

0  242  0  4(1+H)(z-1) 

-24Z  0  a  0 

-24zsin(z+l  )u  24ZCos(z-«-l  )6)  -4(  1+H)  (z-1  )sin(z-l  )£i)  4(  1+H)  (z-1  )cos(z-l  ju 

-24Zcos(z+1  )w  -24zsin(z+l  )u  Acos(z-l)w  AsinCz-Dw 


where 


(A. 14) 


A  =  (24+A)  +  Az  -  H(z-1)(24+A) 


For  the  existence  of  non-trivial  solutions  it  is  necessary  and  sufficient 


det(Z)  =  -8z^4^{(A-4(l+H)(z-l))^sin^u  -  (A+4(H+1 ) (z-1 ) )^sin^zw>  =  0. 


Hence  we  obtain  the  equation 


sin  zw 


rA-4(H-»l)(z- 

Ia+4(H+1 ) (z- 


(z-l)V 

r^nj 


sin  u 


By  (A. 5)  and  (A. 14)  we  have 


(A. 15) 


A-4(H-»-1)(z-1) 


which  implies  (A. 12)  Immediately. 


For  z  »  1,  by  (A. 13)  and  (A. 7),  B  satisfies  ZB  =  0  with 
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0 


2m 


0 


0 


-2m  0  2(m+^)  0 

z  = 

-2MSin2ci>  2MCOs2(d  0  0 

-2mcos2u  -2Msin2(i>  2{n+\)  0 

T 

then  det(£)  =  0,  aoid  (0,0,0, 1)  is  the  corresponding  solution  of  EB  =  0. 

^  T 

Accordingly  u  =  (0,1)  ,  which  represents  a  rotation  around  the  origin. 

For  z  =  -1  by  (A. 13)  and  (A.7).  B  satisfies  ZB  =  0  with 

0  2m(1-M)  0  -2m  ' 

-2m( 1-M)sin2u  -2m( l-M)cosww  0  -2m 

Z  = 

0  0  2m  0 

0  0  2m  0 

then  det(Z)  »  0.  Hence  for  z  =  ±1,  the  equation  ZB  has  a  non-trlvial 
solution.  The  null  space  corresponding  to  the  eigenvalue  t)  =  ±i  is  not  our 
interest  here,  and  we  will  not  elaborate  on  it  further,  but  refer  to  [26]. 

For  z  =  0,  by  (A.  5)  and  (A. 8)  B  satisfies  ZB  =  0  with 


0 

0 

1 

0  ■ 

0 

0 

0 

1 

2m(2m+A) 

0 

0 

cosu 

sinu 

0 

0 

-sinw 

cosw 

Obviously,  det(Z)  =0,  rank(Z)  =2,  and  (1,0, 0,0)^  and  (0, 1,0,0)^  are  two 
linearly  independent  solutions  of  ZB  «  0.  Accordingly  the  space  of 
non-trivial  solutions  of  homogeneous  problem  is  spainned  by 

—  T  —  T 

Sj  *  (cos0,  -sine)  ,  »  (sln8,  cose) 
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which  represent  the  treinslation  in  and  Xg  directions. 


Lemma  A.  3.  Let  t>  =  iz  be  an  eigenvalue  of  the  operator  H(t>)  for  the  mixed 
problem  (A. 2)  eind  (A. 3c),  then  z  satisfies  the  equation 


(A. 16) 


.  2  4(l-u)^  z^  ,  2 

sin  zw  =  j - sin  w, 

3-4y  3-4y 


in  the  case  that  (l+-)cos2w  +1=0,  t)  =  ±i  are  the  eigenvalues,  and 
(A. 16)  is  satisfied. 


Proof.  For  the  mixed  boundary  condition  we  have  u|„  -  =  0  3uid  TL_  = 

0— U  0— tt) 

((T  =0.  If  z  *  0,  ±1  by  (A. 6)  auid  (A.  13),  B  satisfies  ZB  =  0 

T0  00  0®<4) 

with 


-2pzsin(z+l )w  24z  cos(z+1)«  -(i(H+l)(z-l)sin(z-l)w  m(H+1 ) (z-1 )cos(z-l )w 
-2zcos(z-l)w  -2pz  sin(z+l)«  Acos(z-l)«  Asin(z-l)u 


where  A  is  given  by  (A. 14).  Then 


detCZ)  =  4p  Hz'^  -  Am(1+H)(z-1)  +  2m(HAz-pz(z-1  )  ( 1+H) ) 

-  2(l+H)z(A-(z-l)(l+H)M)sin^w  -  2(H-1 )zp(A+p(z-l ) ( 1+H) )sin^zu 


For  the  existence  of  non-trivial  B,  det(Z)  =  0.  Hence  we  have 


(A. 17) 


_  4pHz  -  A(l+H)(z-1)  +  2z(HA-p(z-l)(l+H) 
^  2(1-H)z(A+(z-1)(1+H)4) 


-  CH+l)(A-(z-l)(l+H)p)  2 


K+lz-l)(l+l 


sin  (t>  . 
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Due  to  (A.  H)  and  (A.  14)  we  obtain 


(A. 18) 


(h+lHA-M(z-l)(H+l))  _  t: 
(H-l){A+fi(z-l)(H+l))  ■  3-4u 


(A. 19) 


4nHz  -  A(H+l)(z-l)  +  2z{HA-»i(z-l)(H+l))  _  4(l-w)" 
2(H-1)z(A+(z-1)(1+H)m)  3-4u 


(A. 17)  -  (A. 19)  yields  (A. 16). 

For  z  =  1  by  (A. 13)  and  (A. 7)  B  satisfies  ZB  =  0  with 


-2^sin2(«)  2ti  cos2u 


|-2>icos2w  ~2ti  sin2u 


2(4+\) 


then  det(Z!)  =  4fi^[  ( l+^)cos2«  +  IJ.  Similairly  for  z  =  -1 


[-2fi(  1-M)sin2w  2^(  l-M)cos2(i) 


All  \ 

and  det(£)  »  [(1+  -)cos2u  +  1].  Hence  z  *  ±1  are  eigenvalues  if 

(l+^)cos2w  +1=0,  Eind  also  (A.  16)  is  satisfied  in  the  case  that  z  =  ±1  and 
(l+^)cos2w  +1=0.  Therefore  the  eqviation  ZB  =  0  has  a  non-trivial 

f* 

solution  in  that  case. 


For  z  =  0  then,  by  (A. 13)  and  (A.8)  B  satisfies  ZB  =  0  with 
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0  0  -X  ■ 

1  0 
0  2ti(2ti-*-\)cosu>  2ti(2n+\)sinu 

0  -2ji(2ji+A)sirui)  2|i(2fi+X)cosu 

then  det(Z)  =  24(2ji+X)  >  0.  Hence  ri  ~  0  is  not  am  eigenvalue  of  UiT\) .  □ 

Remairk.  The  trainscendental  equations  for  the  Dirichlet,  Neumajin  auid  mixed 
boundary  conditions  were  derived  by  using  the  argument  of  bihammonic  function 
in  [26]  where  displacement  u  =  Cartesian  coordinates  was  used. 

Consequently  the  coefficients  of  the  equation  (A. 1)  ame  not  constauits. 
Therefore  the  simple  argument  of  linear  system  of  ordinary  equations  with 
constant  coefficients  is  not  valid.  Nevertheless  it  shows  that  the 
tramscendental  equations  are  independent  of  choice  of  coordinates. 
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